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Classical  Model  of  Laser-Stimulated  Surface 
Processes:  Energy  Absorption  Profiles  via  the 
Langevin  Equation 

Jui-teng  Lin  and  Thomas  F.  George"*" 

Department  of  Chemistry,  University  of  Rochester 
Rochester,  New  York  14627,  USA 

A  classical  model  system,  consisting  of  a  diatomic  molecule 
chemisorbed  on  a  solid  surface  and  subjected  to  infrared  laser  radia¬ 
tion,  is  presented.  A  set  of  coupled  equations  of  motion  characterized 
by  the  many-body  effects  of  the  surface  atoms  is  reduced  to  the 
Langevin  equations  of  a  two-body  problem  in  transformed  coordinates. 

The  surf ace- induced  damping  factor  and  frequency  red-shift  of  the 
pumped  mode  are  introduced  by  using  a  Wigner-Weisskopf-type  approxi¬ 
mation.  The  asymmetric  forms  of  the  power  absorption  and  the  quan¬ 
tized  cross  section  due  to  the  nonlinear  effects  of  the  anharmonicity 
are  shown.  The  energy  absorption  profiles  (energy  absorbed  vs  time) , 
which  are  universal  for  any  ranges  of  laser  intensity  (i.e.,  as  the 
intensity  changes  the  profile  remains  the  same  provided  the  time 
scale  is  changed  appropriately)  are  plotted  for  different  sets  of 
the  damping  factor  and  the  detuning.  It  is  found  that  much  longer 
time  scales  (microsecond)  are  required  for  low-power  excitations 
than  for  the  high-power  cases.  The  advantage  and  difficulties  of 
the  normal-mode  method  and  the  numerical  method  are  discussed,  and 
a  new  set  of  coupled  equations  in  the  rotating  frame  are  developed. 

The  energy  absorption  profiles  generated  by  a  numerical  method  of 
a  CO/Ni  system  are  shown. 
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I.  Introduction 

The  interaction  of  infrared  radiation  with  species  chemisorbed 
on  a  solid  surface  has  been  extensively  studied  recently.*-13  In 
our  previous  papers,  the  potential  importance  of  laser-stimulated 
surface  phenomena  in  heterogeneous  systems  (migration,  recombination, 
desorption  and  rate  processes)  in  chemical  industry  and  material 
science  was  pointed  out.  The  possibility  of  controlling  surface 
phenomena  by  means  of  low-power  radiation  was  presented  by  a  simple 
kinetic  model  combining  the  laser  rate  equation  and  the  Langmuir 

q 

equation.  Experimental  evidence  has  suggested  that  laser-stimulated 
surface  processes  (LSSP)  may  be  characterized  by  selective  excitations 
and  hence  are  nonthermal  in  nature.1  The  nature  of  LSSP  (selective 
vs  nonselective)  has  been  more  recently  studied  by  a  theoretical  model 
where  the  level  population  dynamics  of  a  multilevel  system  were  quan¬ 
titatively  discussed.  It  was  shown,  in  a  quantum  mechanical  model, 
that  the  time  scales  of  LSSP  in  a  heterogeneous  system  were  much 
longer  than  those  of  a  gas  phase  system,  due  to  the  fact  that  the 

associated  laser-power  of  the  former  system  was  much  lower  than  that 

11  12 

of  the  latter  system.  '  Furthermore,  the  selectivity  of  the  pumped 
mode,  essentially  characterized  by  the  ratio  of  the  pumping  rate  and 
the  energy  relaxation  rate,  was  analyzed  in  terms  of  a  competition 
between  multiphoton  and  multiphonon  processes.  For  high  selectivity, 
a  long  lifetime  (i.e.,  a  weak  damping)  of  the  excited  pumped  mode  is 
required.12'13 

Most  of  our  previous  models  have  been  cast  in  a  quantum  mech¬ 
anical  framework,  using  a  microscopic  Hamiltonian  which  is  based  on 
the  quantizations  of  the  normal  modes.®-13  The  associated  reduced 
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masses,  the  coupling  factors  and  the  driving  force  of  the  quantum 
Hamiltonian  are  transformed  quantities  defined  in  the  normal  coor¬ 
dinates  and  hence  are  not  explicitly  expressed  as  the  functions  of 
the  original  force  and  the  atomic  coordinates  before  the  normal 
transformation.  Hence,  some  of  the  features  of  the  original  physical 
system  are  not  clearly  displayed  in  quantum  mechanical  calculation. 

We  have  thus  undertaken  a  study  of  a  model  system  based  on  a  clas¬ 
sical  Lagrangian  where  the  masses  and  the  driving  forces  are  defined 
in  terms  of  the  real  atomic  coordinates.  From  a  classical  model  we 
may  investigate  more  details  of  the  physical  picture  of  a  real  sys¬ 
tem  (e.g.,  how  the  active  modes  and  the  associated  generalized  force 
may  be  generated  and  singled  out  from  the  bath  modes) ,  which  may  not 
be  readily  available  in  a  quantum  formulation. 

In  the  present  paper,  we  shall  present  a  classical  model  de¬ 
scribed  by  a  Lagrangian  of  a  system  consisting  of  admolecules  chemi¬ 
sorbed  on  a  solid  surface.  The  many-body  problem  due  to  the  inter¬ 
actions  among  the  admolecules  and  the  surface  atoms  will  be  reduced 
to  a  two-body  problem,  and  then  the  absorption  cross  section  of 
the  system  and  the  stored  energies  of  the  pumped  modes  will  be  com¬ 
puted  analytically.  The  nature  of  LSSP  and  the  excitation  time 
scales,  essentially  characterized  by  the  damping  factor  of  the 
pumped  mode,  will  be  discussed  in  light  of  numerical  results  based 
on  a  set  of  classical  equations  of  motion. 

In  Section  II,  a  model  system  described  by  a  classical 
Lagrangian,  including  the  anharmonicities  of  the  admolecule  and  the 
many-body  couplings  of  both  single  phonons  and  multiphonons  is  pre¬ 
sented.  A  set  of  transformed  equations  of  motion  containing  some 
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important  physical  features  is  discussed.  In  Section  III,  the 
Langevin  equation  of  a  reduced  two-body  problem  of  single-phonon 
processes  is  developed,  where  the  many-body  effects  of  the  surface 
atoms  are  replaced  by  a  damping  factor  and  frequency  red-shift  of 
the  pumped  mode.  We  then  obtain  the  power  absorption  and  the  as¬ 
sociated  cross  section  (in  a  quantization  form)  of  the  system  by 
solving  the  Langevin  equation.  The  steady-state  energies  stored 
in'  the  pumped  modes  are  computed. 

The  universal  energy  absorption  profiles  and  the  time  scales 

2  12  2 

of  LSSP  for  arbitrary  laser-powers  (10  W/cm  -  10  W/cm  )  are  shown 
in  Section  IV.  The  discussion  on  the  normal-mode  method  and  a  set 
of  new  coupled  equations  in  the  rotating  frame,  which  enables  us  to 
overcome  some  difficulties  of  the  usual  numerical  methods,  are  pre¬ 
sented  in  Section  V.  Finally,  we  close  the  paper  in  Section  VI 
with  a  summary  and  conclusion  of  the  main  features  of  LSSP. 

II.  Transformed  Equations  of  Motion  of  a  Many-Body  System 

We  consider  a  model  system  (as  shown  in  Fig.  1)  with  a 
diatomic  molecule  (admolecule)  chemisorbed  on  a  solid  surface  and 
subject  to  infrared  laser  radiation.  The  Lagrangian  may  be  written 

as : 

-  <£,  + 1  *i 
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where  nu,  (i»l,2,3,***)  are  the  mass,  the  displacement  and 

the  frequency  of  the  i-th  atom,  respectively,  and  the  interaction 
terms  with  coupling  constants  X..  and  X.^  (between  the  i-th  and 
j-th  atoms)  are  referred  to  as  the  single-phonon  (linear)  coupling 
and  the  p-phonon  (nonlinear)  coupling,  respectively.  Here  the  ad- 
molecule  is  treated  as  an  anharmonic  oscillator  (up  to  the  quartic 
terms)  while  the  surface  atoms  are  treated  harmonically.  The  an- 
harmonicities  e^,  e2i  and  the  couP^in<3  constant  X^2  are  related 
to  the  derivatives  of  the  potential  energies  and,  e.g.,  for  a  Morse 
potential 

4  2 

we  have  -  -  aJDg,  e12  =  7a  De/12,  and  X12  *  2a  Dg.  Similarly, 
the  coupling  constants  and  X,|p^  are  related  to  the  pair  potential 
energy  between  the  i-th  and  j-th  atoms  by 

jjjMM  x»m  t£v^A  _.v  m  __ . 

The  second  term  in  eq. (1)  is  the  interaction  energy  between  the 
admolecule/surface  system  and  the  laser  field, 

R  (*) ®  j.  E-  c^(cot)  CM.  (9^  ;  (5) 

where  q^  is  the  classical  effective  charge  of  the  atoms,  and  8^  is 
the  angle  between  the  linearly  polarized  electric  field  (with 
circular  frequency  iu)  and  the  coordinate  vector  for  the  optical 
active  mode(s)  of  the  system. 

The  equations  of  motion  associated  with  the  general  form  of 
the  Lagrangian  given  in  eq.(l),  in  principle,  may  be  solved  by 
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the  normal  mode  treatment  or  directly  solved  by  computational  method. 
However,  neither  of  the  above  methods  is  tractable  due  to  the  com¬ 
plicated  many-body  effects  of  the  surface  atoms .  For  the  purpose 
of  some  analytical  results,  we  shall  consider  a  less  general  La- 
grangian,  viz.,  there  are  no  explicit  interactions  assumed  among 
the  surface  atoms  and  the  interaction  energy  between  the  adatoms 
is  assumed  to  be  proportional  to  the  product  of  the  displacement 
from  their  equilibrium  positions. Under  these  assumptions  and  the 
orthogonal  transform 


Ci.t> 

(6-t> 

%rQi>  j*3- 

(6-c) 

the  Lagrangian  of  the  model  system  becomes,  in  the  transformed  normal 
mode  coordinates  (i *  1, 2 ,3, • • • ) , 


Hi 

l  (4  'A'?  )  +  A  Q,  4  +  Oj^A.O. 


+  q£I  Af f  ♦  n 

W  1  r>r‘ 


m'44  +  , 


where  the  transformed  frequencies  are  given  by 
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the  transformed  masses  by 

Mt,i  ”  ( ^  AS.-i  *  ^  4 

Ms  -  'r»i  ,  i  »  3  , 

M”  «  AyA2j  fa -**>/>£  , 


C&a) 


r<p.o 

(M) 

G.e) 


the  transformed  new  anharmonicities  by 

w  _  3[f  5)  ^2j,ij  •+  ~)  -f~  Xa  (A?j  A|j  y  Aw  /\2,-  )J 

4-  ij.  * 

V  '  —  ,  ij  +  £jj  Aij,  2,' 

»V  - t? - , 

amd  the  trains  formed  generalized  forces  by 


a.p 


(fA) 


jj lt}  -  FjW  ,  J  ?  3  .  Cf-il 

The  new  coupling  constants  (note  -  these  are  surface-atom  site 
dependent)  are  defined  by 


7  A I)  Aj; 
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> 


Aj  ~  ( Al]  *  A,j  ^ 

A<f  4-  Al/'Ai 

/”  ( ftll! 

\H>>,  W  ,  . 

Aj  Azi  Z  A*j  Ai> 

J  C'fjt  i ) ! 


r/^.4) 

rffto 

(10J) 


In  deriving  eqs . (7) - (10) ,  we  have  neglected  the  high-order  anharmonic 
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terras  (Q^,  n>  4)  and  considered  the  linear  coupling  terms  between 
Q1  and  Q2  while  keeping  the  high-order  couplings  among  the  adatoms 
and  the  surface  atoms.  Moreover,  the  couplings  among  the  surface 
atoms  (Q^Qj ,  i,  j >  2) ,  which  give  rise  to  an  infinite  number  of 
coupled  equations  of  motion,  are  effectively  absorbed  into  the  site- 
dependent  coupling  constants  (A^,  A2j)  and  tiie  fre<3uencY  dispersion 
of  the  surface-phonon  inodes  (will  be  characterized  by  the  phonon 
mode  spectrum) . 

The  corresponding  equations  of  motion  in  the  transformed 
normal  coordinates  are1^ 


A  1 7JL\  . 

—  I — t~)  “  —  0  ;  »  1  2  2  *  *  ■  ‘ 

dr*Q/  ’ 

and  by  substituting  from  eq. (7) , 


tr1 


(17. CL) 

(12.4.) 


Q2  ♦  ft 

Ql  +S2t°>  -  A  +  Irffl <>z +  Z H 


The  above  equations  of  motion  describe  an  admolecule/surface  system 
with  normal  frequencies  ft2  and  subject  to  the  generalized 

forces  f^(t),  f2(t)  and  f ^  (t)  ,  respectively.  The  important  features 
of  these  transformed  equations  of  motion  are  (see  Fig.  2)  : 

(1)  the  transformation  [eq.(6)l  eliminates  linearly  coupling  between 
the  mode  and  the  surface  phonon  modes  [i.e.,  no  X ^ Q j  in  eq.(12.a)] 

(2)  the  mode  strongly  coupled  to  the  mode  (via  the  AQ2  term)  , 

is  almost  singled  out  of  the  low  frequencies  surface  phonon  modes, 
since  there  is  no  single-phonon  coupling  in  eq . ( l 2 . a )  and  the  p- 
phonon  couplings  [via  are  much  weaker  processes  (based  on 

the  concept  of  the  energy-gap  law);13'16  (3)  from  the  expressions  of 
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the  generalized  forces  [eq.(9)l,  ve  note  that  the  transformed 

applied  field  for  the  case  of  Xi j ~X2 j  and  q 2~_ql '  suggests 

that  the  mode  of  the  admolecule  [corresponding  to  an  asymmetric 

mode  -  see  eq.(S.a)]  may  be  selectively  excited  without  significantly 

heating  the  whole  system  by  means  of  a  radiation  field  with  frequency 
12 

weft,.  We  also  note  for  the  situation  of  A,.>>A..  and  q  ~-q 

^ 3  -*0  *  1 

that  nay  be  comparable  to  f^,  and  hence  either  ft^  or  mey  be 
optically  active  depending  on  the  field  frequency  oj-ft^  or  <jj;;ft-. 

III.  Langevin  Equation  of  Reduced  Two-Body  System 

The  difficulties  of  obtaining  exact  or  analytical  solutions 
of  the  coupled  equations  [eq. (12) ]  are  twofold:  (1)  the  anharmonic 
terms  of  the  equations  of  motion,  in  general,  will  lead  to  time- 
dependent  nonlinear  coefficients  in  the  second  order  differential 
equations;  (2)  the  many-body  effects  of  the  surface  phonon  modes, 
characterized  by  the  site-dependent  coupling  factors,  A  ^  ,  A^^ 
etc.,  will  rule  out  any  tractable  results  when  the  dimension  of  j  is 
large.  In  this  section  we  shall  use  an  iterative  scheme  to  investi¬ 
gate  possible  solutions  of  the  coupled  equations  and  in  turn  compute 
the  power  absorption  and  stored  energy  of  the  system.  We  shall 
first  try  to  linearize  the  anharmonic  terms  and  then  treat  the  many- 
body  couplings  by  some  physically  reasonable  approximations.  As 
mentioned  in  the  previous  section,  the  multiphonon  couplings  (char¬ 
acterized  by  Aj^Q?)  are  in  general  much  weaker  than  that  of  the 
single-phonon  (characterized  by  the  linear  term  AjQ^).  This  may  be 
realized  by  the  fact  that  the  p-phonon  coupling  constant  a(^ 

[see  eqs.(4)  and  (10)]  is  a  strongly  decreasing  function  of  the 
multiphonon  order  p.*^  We  shall  now  examine  the  single-phonon  processes 
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and  neglect  the  multiphonon  terms  for  the  purpose  of  a  simple  treat¬ 
ment  of  the  many-body  effects. 

By  employing  the  asymptotic  (or  harmonic  balance)  method  of 
17 

Bogoluikov,  the  nonlinear  coupled  equations  of  motion  [eq. (12) 
with  no  multiphonon  couplings]  may  be  linearized  as  follows18 


CJ*  +  A  o,  =  [aQj  +  f,tf>  1/m,  y  (isa) 

%  f  =  f  AQ'  *  £XiQj  T  f’*>  J/H  ,  (/3.^ 


We  have  introduced  the  effective  frequencies  and  n2  which  are 
approximately  related  to  the  anharmonicities  (k£ ,  K2)  and  the  steady- 
state  amplitudes  of  the  modes  (A^^,  A2)  by 
-  .*  -2. 


where 


-^1,2  -^1,2  K,2  i 

*  *  5K,». 

f’2  12^  SA.i  • 


In  eq.(13.c),  the  field  acting  on  the  Q^-mode  has  been  dropped  since 

the  low-frequency  surface  phonon  modes  (j >  3)<<  w,  which  are  far 

off-resonant]  are  not  infrared  active.  In  eqs.(13.a)  and  (13. b) ,  the 

terms  of  (M  /M2)Q2  and  (M  are  also  neglected  since  ,Ma  for 

<  X2j  and  m1zm2  (referred  to  eq.(8)]. 

We  shall  next  use  a  technique,  which  enabled  us  to  reduce 

a  multi-level  system  to  a  few-level  system  in  our  previous  quantum 

11  12 

mechanical  models,  '  to  reduce  the  many-body  classical  problem 


to  a  few-body  problem.  To  establish  the  iterative  scheme,  let  us 
represent  the  zeroth-order ,  first-order  homogeneous  solutions  of  the 
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coupled  equations  by  and  Q^,  respectively,  and  choose  the  initial 

conditions  to  be:  Q^(0)  =  Aq,  Q2(0)  =  Bq,  Qj(0)=Aj  (j»3),  Qi(0)  =  0 

(for  all  i) .  By  eq.  (13. b)  we  get 


Qp*)=  6#  o*-(n.2£) }  Os) 

and  using  eq.(15)  we  find  the  homogeneous  solution  of  eq.(13.c)  as 


Q.  (±)~  ^  o (/fcj 

Decomposing  (t-t*)  by  eq.(15)  and  its  derivative,  we  have 

(?(*.*>-  Q"'ci)c^n,t')  -  -k-  .  (0) 

i  *  Sli 

Substituting  eq.(17)  into  eq.(16>  and  working  out  the  integrals  by 
approximating  the  upper  limit  from  t  to  infinity,  we  obtain 

~  0  C~Lriit  T  P  Oij-fO  Q<>  (i} 

where  A  .  »  n_  -  fl.,  At  =  +  fl .  ,  and  P  denotes  the  principal  part. 

3  ^  1  1  *  j 
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In  arriving  at  eq. (18) ,  we  have  used  the  relation 


Substituting  eq. (18)  into  eq.  (13. b)  and  using  the  Wigner- 

19 

Weisskopf-type  approximation,  i.e.,  replacing  the  sum  over  the 
phonon  modes  by  the  integral  over  the  associated  phonon  mode  density 
of  states  p(ftj),  we  find  the  first-order  solution  of  eq.(13.b)  to  be 
embodied  by  the  Langevin  equation: 


where  <••>  denotes  the  ensemble  average  over  the  surface  temperature, 


12 


and  <f  >  is  the  surface  fluctuation  force  given  by 


=<E  ^  Aj-  >  . 


H 


-) 


(21) 


Y  and  Sw,  the  damping  factor  and  frequency  shift,  respectively,  are 
related  to  the  coupling  constant  !■ and  the  phonon  density  of  states 


P  by 


,2  _ 


a  d 

The  above  classical  results  teq.(22)]  are  in  exact  agreement  with  our 


(22,  a.) 
(22.6) 


previous  quantum  mechanical  results  where  the  level  broadening  and 

the  level  shift  correspond  to  the  damping  factor  and  the  frequency 

9  12 

shift,  respectively.  '  It  is  important  to  note  that,  in  eq. (22. a) , 

both  the  coupling  constant  A^{ sT2 )  and  the  phonon  mode  density  p  (?T2 ) 

are  evaluated  at  the  frequency  of  the  n2-mode  which  is  coupled  to  the 

2 

surface  phonon  modes.  For  a  Debye  model  spectrum  p(flj)  * 


with  the  cutoff  frequency  ftD,  we  obtain 

r  =  3icA}(4vfcM,i^nJ)  , 


(2).  a) 


■■ 


3  hi  ^ 


■flj)  + 

“A 


] 


(2U) 


We  also  note  that  in  eq.(20)  the  new  frequency  is  red-shifted  to  the 

t—2  \*4  _  — 

lower  frequency  2  "  since  /«*»  is  a  positive  quantity  for  ^2<^d 

[see  eq.(23.b)].  For  instance,  the  frequency  of  isolated  CO, 

^2  =  2145  ern”^-,  may  be  red-shifted  to  1932  cm  ^  and  2069  cm  ^  when  it 
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is  chemisorbed  on  a  nickel  surface  with  on-top  site  and  bridge  site, 

2i 

respectively . 

So  far,  we  have  reduced  the  many-body  problem  to  a  two-body 
problem  described  by  the  equations  of  motion  (13. a)  and  (20) ,  where 
the  surface-induced  damping  factor  and  frequency  shift  of  the  ft2-mode 
are  introduced  through  the  Langevin  equation.  We  shall  now  solve 
the  equations  and  compute  the  power  absorption  and  the  stored  energy 
of  the  pumped  mode(s).  For  the  case  of  low  surface  temperature 
(i.e.,  )  <f  >r  0,  the  general  solutions  of  the  coupled  equations 

S 

(13. a)  and  (20)  are  found  to  be  in  the  form 

<q®>\  *  «;  A, 

]*=  £  ]  (-24-) 

(QP>I  Ai  hi)  , 

where  R^  are  the  roots  of  the  equation 

rx3*<.fr$)Xi+rn*x *  AAf- =«. 

The  general  solutions  are  complicated  due  to  their  transient  parts. 
However,  for  a  sufficient  time,  the  transient  solutions  vanish,  and 
the  motion  of  the  modes  follows  the  frequency  of  the  field  with  the 
steady-state  solutions  given  by 

>4  H  4)2  \  ( 

W 

Air  An  J  C&lCcAH)  J  2 


(2.7) 

X),=  rwA(M4C?,^+/\^)/(MlWa)  ; 

(2?.  A) 

v2l  -  r*3C,(MAVx+*V,)/Ml  ^ 

W.<) 

I  <qi)> 

I  is. 

\  <  Q2tf)> 


where 


a*-  * cv’^  '*'> 

ft-  <*« 

V-(r“C>/[(*kf-^f,  f*M 

Cf,-Jl*-",  <*4 

nr(^i-^f-  w 

and  the  applied  field  coupling  terms  are  [see  eqs.(5)  and  (9)] 


vr»[A^IE1^,)-Aai,E2^)j/{>aH) ; 


(*>•*) 


Combining  eqs.(9)  and  (30)  and  the  derivative  of  eq.{27)/  we  obtain 
the  steady-state  instantaneous  power  absorption  {force  times  velocity) 
of  the  total  system  (admolecule/surface)  via  the  optically  active 
modes  -  A  mode  (with  frequency  25^)  and  B  mode  (with  frequency  J^)  as 


follows : 


<  Ptf)>  =  t  (  $%>  =<&'*■'>>+<  p;jb) >  '  *'•*-) 


where 


<R m>\  I V, 'A, 


<$/»>)  1  V5A„  -wL MJwt) 


Since  the  sin(2ut)  term  vanishes  as  a  result  of  time  averaging  over 
the  period  of  the  field,  we  immediately  obtain  the  steady-state  average 


power  absorption 


(H.&) 
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The  corresponding  classical  absorption  cross  section,  defined  by 
[power  absorption]/ [laser  intensity],  is  given  by  (in  cgs  units) 

.C&li  ,  o ....  \  (ii .  Q. ) 


(r£  =  Ip-  (<?*»>- <?#») 


Upon  quantization  this  yields  the  quantum  absorption  cross  section 
for  a  transition  between  levels  n  and  (n+1)  of  a  quantum  system*0'18 

V  (*■><)  fl 


<r 


oz 


I 


(Hi*) 


L  [A,-i £*(-»*  *)J2  ~rz2 

where  is  the  detuning  ^  -  u  and  e?  is  the  anharmonicity . 

Note  that  the  absorption  cross  section  is  not  a  symmetric  Lorentzian 
due  to  the  nonlinear  effect,  2e*(n+l/2),  and  the  optimum  detuning  for 
maximum  absorption  cross  section  is  laser  intensity  dependent  (see 
Ref.  10  for  more  rigorous  discussions)  . 

Using  the  steady-state  solutions  [eq. (27) ] ,  we  may  easily 
obtain  the  average  stored  energies  (for  steady-state)  in  the  pumped  A 
and  B  modes,  given  by  EA  and  Eg,  respectively,  as  follows: 

7 


\  -£(AW)(A NAt)  , 

£  =  %■(%**)')  (A»->  An)  . 


(if  jZ) 

i'if.c.) 


For  the  case  of  A  =  0  and  <d  =  n2»  the  stored  energies  reduce  to  the 
simple  forms 


(Sid) 

(}*.*) 


This  is  the  situation  of  very  weak  coupling  between  the  A  and  B  inodes 
where  the  B  mode  is  almost  isolated  from  the  A  mode  but  coupled  to 
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the  bath  (surface  phonon)  modes  via  the  deunping  factor  y.  We  note 
that  for  the  weak  deunping  case  (i.e.,  the  pumped  B  mode  has  a  very 
long  lifetime) ,  one  may  selectively  excite  the  B  mode  without  signi¬ 
ficantly  "heating"  the  other  modes.  A  more  quantitative  description 

of  the  dynamical  nature  of  the  selective  and  nonselective  excitations 

11  12 

has  been  recently  presented  for  a  quantum  mechanical  system.  ' 


IV.  Oniversal  Energy  Absorption  Profiles  _  _  _  _ 

In  the  previous  section,  we  have  reduced  the  many-body 
problem  to  a  two-body  problem  where  the  many-body  effects  are  re¬ 
placed  by  including  a  damping  factor  and  a  frequency  shift  in  the 
Langevin  equations .  In  the  two-body  problem,  the  general  solutions 
of  the  coupled  equations  of  motion  are  still  intractable  due  to  the 
non-explicit  forms  of  the  transient  solutions  [eq.(24)J.  Instead  of 
evaluating  the  steady-state  stored  energies  (eq.(34)],  we  shall  now 
investigate  the  time  evolution  of  the  energy  absorption  of  the  pumped 
mode  by  further  reducing  the  two-body  problem  to  a  single-body  prob¬ 
lem.  For  this  purpose,  we  consider  the  situation  where  only  the  B 
mode  is  resonantly  pumped  (u:  ffj)  being  weakly  coupled  to  the  A  mode 
(A<<  fj)  but  strongly  coupled  to  the  surface  phonon  modes  (via  the 
damping  factor  y) .  The  Langevin  equation  [eq.(20)],  neglecting  the 
surface  fluctuating  force  <f  >,  becomes  the  equation  of  motion  of  a 

9 

damped  enharmonic  oscillator  subject  to  a  generalized  force  f2(t)  * 
[H^cos  (wt)  , 

< T(Q2(^f  DfiOpY- \ ,  (3d 


The  complete  solution  of  eq.(35)  is  found  to  be 
<  Qj(£)>  *  A*i 

t  e  [ A,  +■  L  I  ; 


A) 
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with  the  initial  values  Ag  *<Q2(0)>,  BQ  *  [<Q2  (0)>  +T<Q2  (0)>/2  J/tuj »  where 


*i'«  (nt-(r/J  J* 

07.1) 

Ki  -  rAK  /z  , 

cni) 

Ad  -  , 

Ole) 

Z  *  t(rcd)f^ 

cn^ 

A?*  . 

07*) 

The  constants  and  Afil  are  referred  to  as  the  absorptive  and  the 


elastic  amplitudes  because  the  time-averaged  power  absorption  is 
entirely  due  to  the  out-of-phase  displacement  A^  sinfwt)  [which 
leads  to  an  in-phase  velocity  with  respect  to  the  driving  field 


V2  cos (wt) ] .  The  corresponding  stored  energy  in  the  pumped  B  mode  is 


-Yt 

Eq  |  1  +  e  -;e  2 


Ma) 


where  Eg  is  the  steady-state  energy  given  by 


c  —  JL_  C2  r 

Si  l 


cite) 


32M*  ^  >•  Ay  +  ( r/a  f  J 

Here  we  have  introduced  an  effective  electric  field  acting  on  the  B 


mode  [see  eq.(9>],  for  the  classical  effective  charge  e  »q^»-q2( 


P  AijEimfQi)  -  AijE.CMfflQ 

"  v  1  \  1  \ 


(a,;  -  a-)* 


on 


and  the  optimum  detuning 


-  k!(A2  +  A2,)  which  is  laser 

i  aJD  el 


intensity  dependent  [see  eq.(14)] 
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By  using  I (laser  intensity)  • 


E 2**/(BK/C),  the  steady-state  energy  may  be  expressed  in  a  conven¬ 
er! 


tional  form 


r 


eliy _ 

E0=/.5*.O  i  T\J\  a2.  +  ( r/if  J 


v 


(  .1 
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where  the  units  used  are:  E0(eV),  I(W/cm2),  M2(amu),  e (4 . 8 x  10"10esu) , 
and  both  the  detuning  and  damping  factor  are  in  the  units  of  cm-^ . 

For  the  case  of  y»0,  we  note  that  the  solution  of  eq.(35)  is 

*  ^  =  ^  ^  f  J  ,  (41) 

and  the  corresponding  energy  absorption  of  the  pumped  mode  is 


which  is  an  oscillatory  function  since  the  available  energy,  for  the 
isolated  B  mode,  will  be  necessarily  transferred  back  and  forth  be¬ 
tween  the  pumped  mode  and  the  laser  field  (via  absorption  and  stimu- 

22 

lated  emission,  in  "quantum  mechanical"  language) .  Note  that 

2  2  2 
eq. (42)  reduces  to  EA(t)  *  (Ke  /CM2)It  which  is  proportiontal  to  t 

for  the  exact  resonance  case. 


The  energy  absorption  given  by  eq.(38)  is  shown  in  Fig.  3 
for  different  sets  of  the  optimum  detuning  ^opt  and  the  damping  fac¬ 
tor  y.  It  is  important  to  note  that  these  energy  absorption  profiles 

2 

are  universal  for  all  ranges  of  the  laser  intensity _ (1*10  W/cm  - 

12  2  -1 
10  W/cm  )  when  the  associated  time  scales  in  units  of  y  are 


chosen.  From  eqs.(38)  and  (40),  we  may  define  the  time  scales 
(in  units  of  y”1)  by  relating  the  laser  intensities  (1^  and  I2)  and 
the  damping  factors  (y^  and  yj)  by 


4  3) 


For  instance,  curve  CD)  in  Fig.  3  describes  the  time  dependence  of  the 
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2  2 

energy  absorption  of  the  pumped  mode  for  the  low  power  case  1^  =  10  w/cm 

with  the  time  scale  =  1.3xlo“7,  as  well  as  for  the  high  power 
8  2 

case  I2 *  10  W/cm  ,  but  the  time  scale  is  reduced  to  [by  eq.(43)] 

x  10  ■*.  It  is  seen  from  these  universal  energy  absorption 

profiles  that  the  energies  reach  the  steady-state  values  in  the  micro- 

second  time  scale  for  the  low  power  case  (I-  10  W/cm  )  while  being 

8  2 

in  the  nanosecond  region  for  the  high  power  case  (I  -  10  W/cm  )  and 

12  2 

picosecond  for  I  -  10  W/cm  .  This  is  one  of  the  important  features 

of  laser-stimulated  surface  processes ,  where  low-power  radiation  _ 

2 

(I  =*  10-100W/cm  )  may  be  used  to  study  the  adspecies/surface  system 
in  a  much  longer  time  scale  compared  to  that  of  a  gas-phase  system 

11  12 

(e.g./  multiphoton  dissociation  of  polyatomic  molecules  like  SFg) . 

It  should  also  be  noted  that  for  selective  excitations  by  means  of 

low-power  radiation  to  be  possible,  one  requires  not  only  a  long 

lifetime  of  the  pumped  mode  (i.e.,  small  damping  factor)  but  also 

a  well-defined  laser  frequency  such  that  the  optimum  detuning  AQpt 

-3  -4 

and  the  damping  factor  y  both  have  the  small  values  like  10  -  10 

cm  1  [see  eq. (40) ] . 

For  a  comparison  of  the  energy  absorption  profiles  given  by 

the  reduced  single-body  Langevin  equation  [eq.(35?T  and  those  of  a 

set  of  coupled  equations,  we  show  in  Fig.  4  the  numerical  results 

obtained  by  solving  the  classical  equations  of  motion  are  for  a  model 
23 

system  of  CO/Ni .  We  see  that  the  energy  absorption  curve  shown  m 
Fig.  4  is  different  from  curve  (D)  in  Fig.  3  by  the  fluctuation 
(broadening)  feature  of  the  energy  absorption.  In  the  single-body 
problem,  we  obtain  only  the  average  value  of  the  fluctuating  energy, 
whereas  by  directly  solving  a  many-body  problem  we  may  investigate 


20 


in  detail  fast  oscillations  of  the  energy  absorption  profile  caused 

by  the  energy  relaxation  of  the  pumped  mode  and  the  feedback  from 

the  surface.  We  also  show  the  energy  absorption  profile  for  the 

CO/Ni  system  with  ^Qpt  “Y  s 0  in  Fig.  5.  We  see,  except  for  the 

fluctuating  behavior,  the  resonance  excitation  curve  of  eq.(42) 

2 

EA(t)  *  It  ,  as  expected  in  the  single-body  problem. 


V.  Discussion 

In  this  section  we  shall  discuss  some  advantages  and  difficul¬ 
ties  of  the  normal-mode  method  and  present  a  set  of  trans formed 
equations  of  motion  (in  the  rotating-wave  approach) ,  which  are  more 
practical  in  the  classical  trajectory  calculations  for  the  case  of 
low-power  excitation  processes  where  a  very  long  time  scales  are 
involved. 

A.  Normal-Mode  Method 

Consider  a  Lagrangian  in  the  general  form 

j,  =  Z  ( -  v>i  *  R, ) + 1  ^ ,  U4) 

where  ^ and  are  the  kinetic  energy,  the  dissipation  func¬ 
tion  and  the  potential  energy,  respectively: 


i 

J- 


(4£a) 


(4S.  ■&) 
(tfc) 


The  corresponding  equations  of  motion  of  this  dissipative 
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system  are  given  by 
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3<?C  __  9  9  Ej 

Xt 


(U) 


The  normal-mode  method  capitalizes  on  the  fact  that  each  equation 
of  motion  in  eq. (46)  in  the  normal  coordinates  involves  only  a 
single  coordinate  and  all  the  variables  are  completely  separated. 
However,  the  mass  coefficient  nu  ^  and  the  dissipation  factor 
are  in  general  functions  of  the  coordinate  Xj_,  and  hence  a  set  of 
normal  transforms  which  simultaneously  diagonalizes  T^ ^ ,  F^j  and 

is  not  in  general  available.^-5  For  a  physical  system  consisting 
of  N  adatoms  chemisorbed  on  a  solid  surface,  there  are  3N  "frustrated 
surface  normal  modes  so  that  the  normal-mode  method  essentially 
involves  the  diagonalization  of  a  3N x  3N  matrix,  which  is  in  general 
not  available  particularly  for  the  case  where  the  mass  coefficient 
nu  j  is  coordinate-dependent  and  when  N  is  a  big  number.  For  a 
tractable  model,  we  shall  seek  a  method  in  which  the  dimension  of 
a  related  matrix  may  be  significantly  reduced.  Consider  the  case 
where  the  total  system  (subject  to  a  laser  field)  may  be  decomposed 
into  two  (or  more)  almost  isolated  subsystems,  then  each  subsystem 
involves  the  diagonalization  of  a  small  matrix  which  is  possible  by 

a  set  of  normal  transform.  Such  subsystems  may  be,  in  principle, 

15,(6 

physically  possible  when  there  is  a  big  energy  gap  between  them. 

This  decomposition  may  also  be  mathematically  possible  by  a  set  of 
partially  orthogonal  transforms,  e.g.,  in  Section  I  where  we  used 
the  transformed  coordinates  in  eq.(6)  to  decompose  the  total  system 
into  the  mode  and  a  subsystem  (&2  node  plus  bath  modes)  with  the 
mode  (or/fcnd  the  fi2  node)  able  to  be  selectively  excited. 

For  a  simple  example,  let  us  consider  a  system  which  is  de¬ 
composed  into  two  subsystems,  where  subsystem-2  (with  nor-  1  mo*  is 
referred  to  as  the  inactive  bath  modes)  is  weakly  coupled  to  the 
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subsystem- 1  (with  3N-&  normal  modes  referred  to  as  the  pumped 
mode(s)},  then  the  average  power  absorption  of  the  system  is  given 
by  _  r  .  2 

Pw  =  1 


J2i 


:'4T) 


A.*i 


( ax-uf^(T^)  , 

where  uk  is  a  normal  frequency  of  the  pumped  subsystem  with  the 
associated  damping  factor  and  f ^  is  a  transformed  generalized 
forces  related  to  the  generalized  force  in  the  original  coordinates 
^  (for  driving  forces  given  by  (t)  =  cos  cut)  by 

a 

where  are  the  elements  of  the  transformation  matrix  which  simul- 
tane-usly  diagonalizes  the  kinetic  energy,  dissipation  function  and 
potential  energy  of  subsystem-1,  but  does  not  necessarily  diagonalize 
those  of  the  subsystem-2. 

B.  Numerical  Method  in  the  Rotating  Frame 

Consider  a  model  system  consisting  of  adspecies  (adatoms-1 
and  2)  chemisorbed  on  a  solid  surface  and  subject  to  an  external 
field  [V^  sin(wt)],  described  by  a  set  of  coupled  equations  (for 
only  the  nearest  interactions  are  considered)  as 


%  =  -  ^l7—  -nVi,  ■+ 


*  *  's  ^ 


- 


(4 4) 

Ift) 

i 


where  (i=l ,  2 , 3,  •  •  - )  is  the  coordinate  of  the  i-th  atom  for  the 
longitudinal  motion,  and  the  damping  terms  are  included  to 

take  into  account  the  lateral  interactions  between  the  atoms  of  one 
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row  with  those  of  another.  The  damping  factors  Yj_  { i =1 , 2 )  of  the 
adatoms  [which  may  be  expressed  in  the  form  of  eq.(22.a)]  simulate 
the  surface  effects  of  the  solid  crystal,  and  the  damping  factors 
Yj(j>3)  of  the  solid  atoms  simulate  the  effect  of  the  bulk  of  the 
lattice  in  presenting  the  free  translational  motion  of  the  one¬ 
dimensional  linear  chain.  This  is  the  significant  difference  of 
eqs . (49) - (51)  from  that  of  the  usual  one-dimensional  chain  model, 
where  the  latter  loses  all  the  many-body  surface  effects  of  the 
adatoms  and  the  many-body  bulk  effects  of  the  solid  atoms. 

As  discussed  in  Section  IV,  the  time  scales  of  the  energy 

absorption  profiles  (in  units  of  y  ^)  for  the  low-power  excitations 
2  2 

(I  -  10  W/cm  )  are  in  the  ranges  of  microseconds,  which  are  much 

-14 

longer  than  the  oscillation  cycle  of  the  field  (w  -  10  second)  . 

This  causes  the  difficulty  in  obtaining  the  absorption  profiles  by 
a  computational  method  which  solves  the  above  coupled  equations 
directly.  Furthermore,  the  energy  absorption  of  the  pumped  mode 

2 

is  very  sensitive  to  the  amount  of  detuning  (A  =  tOg-w)  .  For  I  -  10 
W/cm2,  we  require  that  AQpt,  y  *  1Q-3_  1(J-4  ^-1  which  also  causes 

difficulty  in  tuning  the  field  frequency  to  obtain  optimum  excitations 

8  2 

Note  that  for  the  cases  of  high-power  excitations  (I>  10  'W/cm  ), 
these  difficulties  would  not  be  encountered  and  hence  we  may  obtain 
the  absorption  profiles  numerically  from  eqs . (49) - (51)  [Figs.  4  and 
52. 

In  order  to  overcome  the  difficulties  described  above-  and 
obtain  the  energy  absorption  profiles  by  a  numerical  method,  par¬ 
ticularly  for  the  low-power  excitation  processes,  we  now  present 
a  method  which  relies  on  the  rotating  wave  approximation  (KWA) , 
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IQ  y  y 

a  well-known  technique  for  quantum  mechanical  systems.  '  By 
using  the  rotating  frame 

where  Wq  is  the  frequency  of  the  optically  active  mode  of  the  ad- 

species/surface  system,  we  shall  consider  the  near-resonant 

excitation  with  the  detuning  A  =  oJq-oj  z  0  (for  the  harmonic  model) 

*  2 

or  -  oisO  [for  the  anharmonic  model  -  see  eq.  (  14)  ]  . 

The  coupled  equations  of  motion  (49) -(51)  become 

-  -  .2-.  /  „*  .  •>.  „.  \  .  "W 


V « '  » + 1 


(&) 


M) 


e'*at>  i9i>  ISS) 


where  A  =  uiQ  -  ui  is  the  detuning.  In  obtaining  eqs .  (53) - (55)  we  have 
used  the  RWA,  that  is,  we  have  neglected  the  fast  oscillating  terms 
exp  [±  (oJq+w)  t]  .  The  inportant  features  of  the  new  equations  of  motion 
are:  (1)  due  to  the  complex  coefficients,  the  j  coupled  equations 

are  in  fact  2j  real  equations  (equivalent  to  4j  first-order  differ¬ 
ential  equations  which  must  be  numerically  solved) ;  (2)  the  time- 
dependent  field  with  the  very  fast  sinusoidal  function  sin(wt)  is 
eliminated  in  the  rotating  frame  in  the  RWA  and  the  coupled  equations 
are  characterized  by  the  detuning  A  which  leads  to  a  much  slower 
oscillating  function  exp(-iAt)  for  near  resonance. 

Therefore,  by  using  the  new  set  of  equations  of  motion,  we 
are  able  to  compute  the  energy  absorption  by  solving  the  coupled 
equations  numerically.  It  is  important  to  note  that  the  new  coupled 


equations  (53) -(55)  may  be  solved  by  using,  for  example,  the  Runge- 
Kutta  method  with  a  much  longer  time  step  (s  10  sec)  than  for  the 
original  coupled  equations  (49) -(51)  (s  10_16sec) .  The  applications 
of  the  new  coupled  equations  on  some  real  adspecies/surface  system, 
e.g .,  co/Ni  and  H/Pt,  are  in  progress. 

VI.  Summary  of  the  Main  Features  of  LSSP 

We  conclude  the  main  feature  of  LSSP  and  summarize  the 
specific  results  obtained  in  this  classical  model  as  follows: 

(a)  By  a  set  of  orthogonal  transforms,  we  are  able  to  generate 
two  normal  inodes  where  the  high  frequency  asymmetric  (ft^)  mode  is 
uncoupled  to  the  bath  mode  while  being  coupled  to  the  low  frequency 
symmetric  (Bj)  mode.  From  the  transformed  frequencies  teq.(_8.a)] 
and  the  generalized  forces  [eq.(S.a)],  we  are  able  to  see  some  selec¬ 
tive  nature  of  the  system. 

(b)  The  nonlinear  coupled  equation  due  to  the  anharmoni cities 
of  the  potential  energies  may  be  linearlized  by  the  asymptotic  (or 
harmonic  balance)  method  where  the  amplitude-dependent  frequencies 
are  presented  in  eq.(14). 

(c)  By  the  iterative  procedure,  the  many-body  effects  of  the 

surface  atoms  are  absorbed  into  the  Langevin  equation,  where  a 

damping  factor  and  frequency  red-shift  are  introduced  [eq.(22)J. 

These  classical  results  are  in  agreement  with  those  of  the  quantum 

mechanical  calculation  presented  in  our  earlier  work,  where  the 

9-12 

Markovian  statistics  and  Wigner-Weisskopf  theory  were  used. 

(d)  The  power  absorption  and  the  quantized  cross  section  are 


obtained  by  solving  the  coupled  Langevin  equations  and  the  asymmetric 
behavior  is  shown  by  the  nonlinear  effects  of  the  anharmoni cities 
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[eqs . ( 32)  and  (33) ] . 

(e)  The  energy  absorption  profiles  of  the  pumped  mode,  which 

2  12  2 

are  universal  for  any  laser-power  (ranging  from  lOW/cm  to  10  W/cm  )  , 
are  shown  for  different  sets  of  the  damping  factor  and  the  detuning 
[Fig.  3] .  The  long  time  scales  of  LSSP  (in  the  range  of  microseconds) 
are  discussed  for  the  case  of  low-power  excitations  [eq.(43)]. 

(f)  The  average  power  absorption  of  a  system,  which  may  be 
decomposed  into  almost  isolated  subsystems,  is  computed  by  the 
normal-mode  method  [eq.(47)l. 

(g)  Finally,  the  difficulties  in  numerical  methods  of  solving 

a  set  of  coupled  equations  are  pointed  out,  and  a  set  of  new  coupled 
equations,  which  enable  one  to  numerically  generate  the  energy 
absorption  curves  of  low-power  LSSP,  are  presented  in  the  rotating 
frame  (eqs .  (53) -  (55)  1  .  - 
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Figure  Captions 

Fig.l.  Diatomic  molecule  chemisorbed  on  a  solid  surface.  The 

coupling  constants  between  the  adatom-1  and  -2  and  among  the 
adatoms  and  the  surface  atoms  are  given  by  X12,  ^ij  and 
X2j,  respectively. 


Fig. 2.  Schematic  diagram  of  the  density  of  states  for  the  transformed 

normal  frequencies  Q2  and  given  by  eq.(8).  The  coupling 
factors  among  different  modes  are  denoted  by  A  and  \j  [given 
by  eq.(10)].  Note  that  the  mode  is  uncoupled  to  fij  modes 
for  the  single-phonon  processes . 


Fig. 3.  The  universal  energy  absorption  profiles  of  the  pumped  mode  for 

different  sets  of  the  optimum  detuning  AQpt  and  the  damping 

-1  -3 

factor  y  (in  units  of  cm  )  :  curve  A  -  AQpt=»0,  y»2xl0  ; 

-4  -  3 

curve  B  -  Aopt=0.5y  *  5x10  ;  curve  C  -  Aopt*y=10  • 

curve  D  -  Aopt3S2 .5y»2 .5x10  4?  curve  E  -  Aopt*5y=2xl0  5; 

—  e  —4 

curve  F  -  A  =8y=10  ;  curve  G  =  A__ . *10y»10  ;  for  low- 

opt  opt 

2 

power  laser  l=100W/cm  .  Note  that  the  time  scales  are  shown 
in  units  of  y”*. 


Fig. 4.  Energy  absorption  of  the  pumped  CO  mode  in  the  .CO/Ni  system  for 
laser-power  I=10^^W/cm^  with  AQpt=0 .2y*20cm 

Fig. 5.  Energy  absorption  of  the  same  system  as  in  Fig. 4,  but  with 


4opt- 1-  0  • 


ENERGY  ABSORPTIONCeV ) 


0  2  *  S  8  10  12  14  16 


TIME  tCIN  UNITS  OF  T"0 
f;3 . 3  (A) 


TIME  t  C PICOSECOND!) 


TIME  t  CPICOSECONCO 


SP472-3/A1 


472 :GAN:716:ddc 
78u472-608 


TECHNICAL  REPORT  DISTRIBUTIOM  LIST,  GEN 


No. 

Copies 

Office  of  Naval  Research 

Attn:  Code  472 

800  North  Quincy  Street 

Arlington,  Virginia  22217  2 

ONR  Branch  Office 
Attn:  Dr.  George  Sandoz 
536  S.  Clark  Street 

Chicago,  Illinois  60605  1 

OKR  Area  Office 
Attn:  Scientific  Dept. 

715  Broadway 

New  York,  New  York  10003  1 

ONR  Western  Regional  Office 
1030  East  Green  Street 

Pasadena,  California  91106  1 

ONR  Eastern/Central  Regional  Office 
Attn:  Dr.  L.  H.  Peebles 

Building  114,  Section  D 
666  Sunnier  Street 

Boston,  Massachusetts  02210  1 

Director,  Naval  Research  Laboratory 
Attn:  Code  6100 

Washington,  D.C.  20390  1 

The  Assistant  Secretary 
of  the  Navy  (RE&S) 

Department  of  the  Navy 
Room  4E736,  Pentagon 

Washington,  D.C.  20350  1 

Commander,  Naval  Air  Systems  Command 
Attn:  Code  310C  (H.  Rosenwasser) 

Department  of  the  Navy 

Washington,  D.C.  20360  1 

4 

Defense  Technical  Information  Center 
Building  5,  Cameron  Station 
Alexandria,  Virginia  22314  12 


No. 

Copies 

U-S.  Army  Research  Office 
Attn:  CRD-AA-IP 
P.0.  Box  12211 

Research  Triangle  Park,  N.C.  27709  1 

Naval  Ocean  Systems  Center 

Attn:  Mr.  Joe  McCartney 

San  Diego,  California  92152  1 

Naval  Weapons  Center 
Attn:  Dr.  A.  B.  Amster, 

Chemistry  Division 

China  Lake,  California  93555  1 

Naval  Civil  Engineering  Laboratory 
Attn:  Dr.  R.  W.  Drisko 

Port  Hueneme,  California  93401  1 

Department  of  Physics  &  Chemistry 

Naval  Postgraduate  School 

Monterey,  California  93940  1 

Dr.  A.  L.  Slafkosky 
Scientific  Advisor 
Commandant  of  the  Marine  Corps 
(Code  RD-1) 

Washington,  D.C.  20380  1 

Office  of  Naval  Research 
Attn:  Dr.  Richard  S.  Miller 
800  N.  Quincy  Street 

Arlington,  Virginia  22217  1 

Naval  Ship  Research  and  Development 
Center 

Attn:  Dr.  G.  Bosmajian,  Applied 
Chemistry  Division 

Annapolis,  Maryland  21401  1 

Naval  Ocean  Systems  Center 
Attn:  Dr.  S.  Yamamoto,  Marine 
Sciences  Division 

San  Diego,  California  91232  1 


Dr.  Fred  Saalfeld 
Chemistry  Division,  Code  6100 
Naval  Research  Laboratory 
Washington,  D.C.  20375 


Mr.  John  Boyle 

Materials  Branch 

Naval  Ship  Engineering  Center 

Philadelphia,  Pennsylvania  19112 


1 


1 


SP472-3/A3 


472 :CAN: 716 : ddc 
78u472-608 


TECHNICAL  REPORT  DISTRIBUTION  LIST,  GEN 
No. 

Copies 


Dr.  Rudolph  J.  Marcus 
Office  of  Naval  Research 
Scientific  Liaison  Group 
American  Embassy 
APO  San  Francisco  96503 

Mr.  James  Kelley 
DTNSRDC  Code  2803 
Annapolis,  Maryland  21402 


Dr.  David  L.  Nelson 
Chemistry  Program 
Office  of  Naval  Research 
800  North  Quincy  Street 
Arlington,  Virginia  22217 


SP472-3/A17 


472 : GAN: 716 : ddc 
78u4 72-603 


TECHNICAL  REPORT  DISTRIBUTION  LIST,  056 


No. 

Copies 

Dr.  D.  A.  Vroom 
IRT 

P.O.  Box  80817 

San  Diego,  California  92138  1 

Dr.  G.  A.  Somorjai 

Department  of  Chemistry 

University  of  California 

Berkeley,  California  94720  1 

Dr.  L.  N.  Jarvis 
Surface  Chemistry  Division 
4555  Overlook.  Avenue,  S.W. 

Washington,  D.C.  20375  1 

Dr.  J.  B.  Hudson 

Materials  Division 

Rensselaer  Polytechnic  Institute 

Troy,  New  York  12181  '  1 

Dr.  John  T.  Yates 
Surface  Chemistry  Section 
National  Bureau  of  Standards 
Department  of  Commerce 

Washington,  D.C.  20234  1 

Dr.  Theodore  E.  Madey 

Surface  Chemistry  Section 

Department  of  Commerce 

National  Bureau  of  Standards 

Washington,  D.C.  20234  1 

Dr.  J.  M.  White 
Department  of  Chemistry 
University  of  Texas 

Austin,  Texas  78712  1 

Dr.  Keith  H.  Johnson 

Department  of  Metallurgy  and  Materials 
Science 

Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139  1 

Dr.  J.  E.  Demuth 
IBM  Corportion 

Thomas  J.  Watson  Research  Center 
P.0.  Box  218 

Yorktown  Heights,  New  York  10598 


No. 

Copies 

Dr.  C.  P.  Flynn 
Department  of  Physics 
University  of  Illinois 

Urbana,  Illinois  61801  1 

Dr-  W.  Kohn 
Department  of  Physics 
University  of  California 
(San  Diego) 

LaJolla,  California  92037  1 

Dr.  R.  L.  Park 
Director,  Center  of 
Materials  Research 
University  of  Maryland 

College  Park,  Maryland  20742  1 

Dr.  W.  T.  Peria 
Electrical  Engineering 
Department 

University  of  Minnesota 

Minneapolis,  Minnesota  5.5455  1 

Dr.  Narkis  Tzoar 

City  University  of  New  York 

Convent  Avenue  at  138th  Street 

New  York,  New  York  10031  1 

Dr.  Chia-wei  Woo 

Department  of  Physics 

Northwestern  University 

Evanston,  Illinois  60201  1 

Dr.  D.  C.  Mattis 
Polytechnic  Institute  of 
New  York 
333  Jay  Street 

Brooklyn,  New  York  11201  1 

Dr.  Robert  M.  Rexter 

Department  of  Chemistry 

University  of  Minnesota 

Minneapolis,  Minnesota  55455  1 

Dr.  R.  P.  Van  Duyne 
Chemistry  Department 
Northwestern  University 
Evanston,  Illinois  60201 


1 


1 


SP472-3/A19 


472 :GAK : 716 : ddc 
78u472-608 


TECHNICAL  REPORT  DISTRIBUTION'  LIST,  056 


No. 

Copies 

Dr.  M.  G.  Lagally 
Department  of  Metallurgical 
and  Mining  Engineering 
University  of  Wisconsin 
Madison,  Wisconsin  53706  1 

Dr.  Robert  Gomer 
Department  of  Chemistry 
James  Franck  Institute 
5640  Ellis  Avenue 

Chicago,  Illinois  60637  1 

Dr.  R.  G.  Wallis 

Department  of  Physics 

University  of  California,  Irvine 

Irvine,  California  92664  1 

Dr.  D.  Ramaker 

Chemistry  Department 

George  Washington  University 

Washington,  D.C.  20052  1 

Dr.  P.  Hansma 
Chemistry  Department 
University  of  California, 

Santa  Barbara 

Santa  Barbara,  California  93106  1 

Dr.  P.  Hendra 
Chemistry  Department 
Southhampton  University 

England  S09  SlM  l 

Professor  P.  Skell 
Chemistry  Department 
Pennsylvania  State  University 
University  Park,  Pennsylvania  16802  1 

Dr.  J.  C.  Hemminger 

Chemistry  Department 

University  of  California,  Irvine 

Irvine,  California  92717  1 

Dr.  Martin  Fleischmann 
Department  of  Chemistry 
Southampton  University 
Southampton  S09  5NH 

Hampshire,  England  1 


No. 

Copies 

Dr.  J.  Osteryoung 
Chemistry  Department 
SUNY,  Buffalo 

Buffalo,  New  York  14214  1 

Dr.  G.  Rubloff 
I.B.M. 

Thomas  J.  Watson  Research  Center 
P.  0.  Box  218 


Yorktovn  Heights,  New  York  10598  1 

Dr.  J.  A.  Gardner 

Department  of  Physics 

Oregon  State  University 

Corvallis,  Oregon  97331  1 

Dr.  G.  D.  Stein 

Mechanical  Engineering  Department 

Northwestern  University 

Evanston,  Illinois  60201  1 

Dr.  K.  G.  Spears 

Chemistry  Department 

Northwestern  University 

Evanston,  Illinois  60201  1 

Dr.  R.  W.  Plummer 
University  of  Pennsylvania 
Department  of  Physics 

Philadelphia,  Pennsylvania  19104  1 

Dr.  E.  Yeager 

Department  of  Chemistry 

Case  Western  Reserve  University 

Cleveland,  Ohio  41106  2 

Professor  George  H.  Morrison 

Cornell  University 

Department  of  Chemistry 

Ithaca,  New  York  14853  1 


Professor  N.  Winograd 
Pennsylvania  State  University 
Chemistry  Department 

University  Park,  Pennsylvania  16802  1 

Professor  Thomas  F.  George 

Department  of  Chemistry 

University  of  Rochester 

Rochester,  New  York  14627  1 


2 


